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Abstract 

The scattering of atoms and rotons at the free surface of superfluid 4 He is 
studied in the framework of linearised time dependent mean field theory. The 
phenomenological Orsay- Trento density functional is used to solve numeri- 
cally the equations of motion for the elementary excitations in presence of a 
free surface and to calculate the flux of rotons and atoms in the reflection, 
condensation, and evaporation processes. The probability associated with 
each process is evaluated as a function of energy, for incident angles such 
that only rotons and atoms are involved in the scattering (phonon forbidden 
region). The evaporation probability for R + rotons (positive group velocity) 
is predicted to increase quite rapidly from zero, near the roton minimum, to 
1 as the energy increases. Conversely the evaporation from R~ rotons (neg- 
ative group velocity) remains smaller than 0.25 for all energies. Close to the 
energy of the roton minimum, A, the mode-change process R + <-> R~ is the 
dominant one. The consistency of the results with general properties of the 
scattering matrix, such as unitarity and time reversal, is explicitly discussed. 
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The condensation of atoms into bulk excitations is also investigated. The 
condensation probability is almost 1 at high energy in agreement with ex- 
periments, but it lowers significantly when the energy approaches the roton 

minimum in the phonon forbidden region. 
PACS number: 67.40.-w, 67.40.Db 
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I. INTRODUCTION 



One of the unique features of superfluid helium at low temperature is that elementary 
excitations like rotons and high energy phonons have a very long mean free path, long enough 
to propagate ballistically in the liquid on macroscopic distances. Furthermore their energy 
is comparable with the energy required to eject atoms from liquid to vacuum. Thus rotons 
and phonons impinging on the free surface of the liquid can evaporate atoms by one-to- 
one quantum processes. One can also produce collimated beams of elementary excitations 
and collect the evaporated atoms above the surface in order to extract information on the 
properties of the superfluid. Similarly one can investigate the opposite process of atoms 
which condense in the liquid by producing elementary excitations. Several experiments 
have been carried out in the last decades fllHlOf to explore the phenomenon of quantum 
evaporation and condensation (see also Ref. [0 for a recent review). From the theoretical 
viewpoint |l2|- |T9fl these processes are very appealing, being clean examples of scattering of 
excitations in a highly correlated many-body quantum system. A good theory of quantum 
evaporation and condensation is however a difficult task and a clear understanding of the 
fundamental mechanism underlying these phenomena has still to come, even though relevant 
steps in this direction have been accomplished in the last years. 

In this paper we present a time dependent mean field theory, based on the density 
functional formalism, which allows one to calculate the scattering matrix elements associated 
with the scattering of rotons, phonons and atoms at the free surface. Preliminary results of 
this theory have been recently published in Ref. [19|]. Here we discuss the theory in more 
detail and we present results in a wider range of energy and wave vectors. 

A discussion about quantum evaporation and condensation requires first a detailed 
knowledge of the spectrum of elementary excitations of superfluid 4 He. A schematic picture 
is given in Fig. [1]. The liquid- vacuum interface is supposed to be in the xy-pl&ne. The system 
is translational invariant in those directions, so that the parallel wave vector q x is a conserved 
quantity (we choose q y = without any loss of generality). The minimum energy needed 
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to produce free atom states outside the liquid at zero temperature is the chemical potential 
= 7.15 K. The threshold for evaporation for given values of q x is then + h 2 q^/2m 
(dashed line). Phonons, rotons and atoms can propagate at different angles with respect 
to the surface, that is with different q z for a given q x . The dispersion law for phonons and 
rotons propagating in the bulk liquid and in the direction parallel to the surface (q z = 0) is 
plotted as a function of q x . Excitations having q z ^ fill different regions of the spectrum 
in Fig. m with a continuum of states. For instance each point in region V corresponds to 



a roton R + (positive group velocity) propagating with q = \Jqi + ql such that fiuj(q) is 
the R + roton dispersion. Since q x < q, roton states with q z ^ fill the region at the left 
of the roton branch plotted in the figure. Similarly in region IV also R~ rotons (negative 
group velocity) can propagate; each point in that region can be a combination of R + and 
R~ rotons propagating at different angles. Above the threshold for evaporation, in region 
III, atom states outside the liquid are possible too. Similarly in region II one has phonons, 
rotons as well as atoms; in region I, below the energy of the roton minimum, A, there are 
only phonons and atoms; finally in region VI, above the maxon energy, there are rotons and 
atoms. The free surface acts as a scattering region for all these excitations. 

In Refs. [0,0] the effect of the surface on the propagation of the elementary excitations 



has been studied within a local density approximation. The dispersion law of the excitations 
was assumed to be the same as for a uniform liquid of density equal to the local density. 
Maris Jl7j used a simple interpolation between the phonon-roton dispersion law in bulk liquid 



and the dispersion law of the free atoms. He studied the behavior of classical trajectories 
for quasiparticles along the z-direction assuming the conservation of energy and parallel 
momentum. From that analysis one finds some interesting constraints on the structure of 
the classical trajectories crossing the free surface. For instance only phonons and rotons 
above the maxon energy (about 14 K) are found to evaporate atoms. Vice-versa, the theory 
predicts no evaporation from rotons with energy smaller than the maxon energy, because 
of the occurrence of a barrier at the interface. The experimental evidence |]TT| for quantum 
evaporation induced by rotons even below the maxon energy is consequently an important 



indication of the crucial role played by quantum effects. Quantum effects were partially 
included in the calculations of Ref. |jTB[ by means of perturbation theory based on WKB 



states; the density profile in that case was approximated by a suitable analytic function 
and the dispersion law for the excited states, in local density approximation, was calculated 
using Beliaev's theory with an effective interaction. 

For an accurate treatment of quantum effects one has to go beyond the local density 
approximation and the semiclassical treatment. To accomplish the task we proceed in two 
distinct steps: 

• We first evaluate the ground state and the excited states of the inhomogeneous system 
within a self-consistent theory sufficiently accurate to reproduce the main features of 
the spectrum in Fig. [I], including the structure and the excitations of the free surface. 

• We identify the current carried by the elementary excitations and use the numerical 
solutions of the equations of motion to calculate the asymptotic flux of excitations 
associated with a given scattering process. While this identification is trivial in the 
description of free atoms in vacuum, it becomes less obvious in the liquid where many- 
body effects, present in this highly correlated system, must be properly included. 
Taking the ratios of incoming and outgoing fluxes one finally gets the evaporation and 
condensation probabilities, which are the final goal of the theory. 

For the first step we use a time dependent mean field theory based on the least ac- 
tion principle applied to an energy functional JdrTC[ty*, ty]. For the form of the functional 
we follow the recent proposal of Ref. [2(J. We consider linear variation of ty, in the form 



(tyo + Sty, and expand Sty in plane waves along the parallel direction. We restrict our- 
selves to variations of the energy linear in Sty. For this reason the theory describes only 
one-to-one processes (evaporation of one atom by one roton, creation of one roton by con- 
densing one atom, etc.) and not processes involving more than two excitations (for example 
multi-phonon or multi-ripplon production in the atom condensation process). It is worth 
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mentioning that some experiments by Wyatt and co-workers P,|6|,|TT| support the idea that 
the evaporation of atoms from bulk excitations is essentially a one-to-one process, so that 
the use of a linear theory seems justified in this case. Conversely measurements of atom 
condensation []3],|9| point out an important role played by non linear effects associated with 
the creation of ripplons. This asymmetry between evaporation and condensation is still 
basically unexplained. Quantitative predictions for the scattering rates in linear theory are 
essential in order to better understand the origin of this and other discrepancies. 

As concerns the second step we derive the definition of flux (or current density) associated 
with the propagation of the elementary excitations, in terms of the linearised solutions of 
the equations of motion. We show that such a current density obeys the proper equation of 
continuity. We will devote detailed discussions about the meaning of the flux of excitations 
and about the properties of the scattering matrix elements, which are calculated starting 
from the knowledge of the fluxes involved in the scattering process. 

The theory can be applied to all the relevant regions of the spectrum in Fig. |1]. In 
the present work we present results for processes involving rotons and atoms in the region 
III. We chose this region for two basic reasons: first, phonons are not allowed there and 
thus the analysis is simpler; second, the arguments based on classical trajectories predict no 
evaporation in this region and this means that quantum effects should play a crucial role 
in the evaporation process. The analysis of the remaining parts of the spectrum will be the 
object of a future work. 

The paper is organised as follows: in Sec. we present the theoretical formalism of the 
density functional and the equations of motion. In Sec. |T| we introduce the concept of 
current of excitations. In Sec. |IV| we discuss general properties of the scattering matrix. In 
Sec. [V] we present the results for the probabilities of evaporation and condensation. Finally 
we will give a brief summary of the main results and discuss some future perspectives. 
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II. TIME DEPENDENT DENSITY FUNCTIONAL THEORY 



The formalism used in the present work was already introduced in Ref. |2U|| . In that 



paper detailed discussions were devoted to the motivations of the theory and to the choice 
of the density functional. Results for static and dynamic properties of superfluid 4 He in 
different geometries were presented. In Refs. |2T|J22|] the dynamics of the free surface and of 
droplets, respectively, were studied in detail within the same scheme. As a consequence we 
refer to the above-mentioned papers for details and we recall here simply the main points. 
The starting assumption is that the energy can be written as a functional of the form 



E = JdrH^^*} (1) 
where the complex function \I/ is written as 



*(r,t) = <f>(r,t)ex V ^-S(v,t)) . (2) 

The real function $ is related to the diagonal one-body density by p = $ 2 . The phase S 
fixes the velocity of the fluid through the relation v = (l/m)VS, where m is the mass of 
the 4 He atoms. In the calculation of the ground state only states with zero velocity must 
be considered, so that the energy is a functional only of the particle density p(r). A natural 
representation is given by 

E = jdvHM = E c [p] + Jdr^-{S7Jp) 2 , (3) 

where the second term on the r.h.s. is a quantum pressure, corresponding to the kinetic 
energy of a Bose gas of non uniform density. The quantity E c [p] is a correlation energy; it 
incorporates the effects of dynamic correlations induced by the interaction. Ground state 
configurations are obtained by minimising the energy of the system with respect to the 
density. This leads to the Hartree-type equation 



(# Q -/i)v7>(r) = , (4) 

where 



H = ~V 2 + U\p,t] (5) 

is an effective Hamiltonian. The quantity U[p, r] = 5E c /5p(r) acts as a mean field, while 
the chemical potential /i is introduced in order to ensure the proper normalisation of the 
density to a fixed number of particles. The dynamics of the system can be studied by using 
the least action principle: 



/t'2 
dt I <!r 



. (6) 



By making variations with respect to ^ or ^* one derives the equations of motion for the 
excited states, in the form of Schrodinger-like equations. The equation for \1/ is 

(1T- M )¥ = ift^* , (7) 

where H = 5E/8^>* is an effective Hamiltonian. We linearise the equation by writing 

tf(r,t) = * (r)+$¥(r,t) (8) 

where \I/o( r ) corresponds to the ground state. The change of wave function 5^> can be written 
in the form 

^(r,t) = /(r)e-^ + /(r)e^, (9) 

where the functions /(r) and g(r) are fixed, together with the frequency lu, by the solution 
of the equations of motion ([I]). The Hamiltonian H then takes the form 

H = H + 5H (10) 

and the equation of motion becomes 

(H - /i)5^ + 5H = ifi^-Sy . (11) 

The term 5H is linear in 5^ and accounts for changes in the Hamiltonian induced by the 
collective motion of the system. Since H depends explicitly on the wave function Equation 
(0) has to be solved using a self-consistent procedure. 
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The above scheme corresponds to a time dependent density functional (TDDF) theory. 
A discussion about equivalent formulations of the same approach (hydrodynamic-like equa- 
tions, Green's functions), as well as about the connection with other many body approaches, 



is given in Refs. p0|-|22|j. Equations of the same type are obtained using Correlated Basis 
Functions in the context of Hypernetted Chain approximation [f23] . Here we stress that the 
theory corresponds to a quantum mechanical treatment of the dynamics and consequently 
accounts for the interference and tunneling phenomena which are expected to play a cru- 
cial role in the evaporation process. Of course, due to linearisation, they do not include 
inelastic processes associated with multi-phonons or multi-ripplons. As already said in the 
Introduction, these effects lie beyond the present theory. 

To proceed further one has to specify both the explicit form of the energy functional 
(H) and the geometry of the system under investigation. The energy functional is the one 



introduced in Ref ||20|| . It consists of two parts, a functional of the particle density only 
and a functional of both the particle density and the current density. It corresponds to a 
generalisation of a previous finite-range functional introduced by Dupont-Roc et al. |]24}| . Its 
explicit form is given in Appendix A. The functional is phenomenological, i.e., it contains 
parameters which are fixed to reproduce known properties of the bulk liquid. In particu- 
lar the equation of state and the static response function in bulk liquid are reproduced by 
construction. The new functional, with respect to the one of Ref. [|4]], contains a current- 
current effective interaction, which accounts for backflow-like correlations and allows one to 
reproduce the experimental phonon-roton spectrum in the uniform liquid. The latter is a 
crucial requirement for a theory whose goal is the prediction of evaporation and condensa- 
tion probabilities, since these quantities depend dramatically on the form of the excitation 
spectrum. 

As concerns the geometry, we assume translational invariance in the x-y plane. The 
density profile of the liquid in the ground state is a function of the orthogonal co-ordinate z. 
The profile of the free surface obtained with the present density functional theory is about 
6 A thick and compares well with the one coming from ab initio calculations as well as with 



the experimental data of Ref. |25j (see Ref. |22| for a discussion). The excited states are 



described by fluctuations of the particle density and velocity associated with the solutions 
(§) of the equations of motion. In the planar geometry the same solutions can be written in 
the form |26| 

6V{r,t) = f{z)e iigxX - ut) + g{z)e- ii(lxX - UJt) . (12) 

With a proper choice for the boundary conditions the functions f(z) and g(z) can be chosen 
real. Once the equations of motion are written explicitly for the unknowns f(z) and g(z), 
they assume the typical form of the equations of the random phase approximation (RPA) 
for bosons. In particular they account for both particle-hole [/(#)] and hole-particle [g{z)\ 
transitions which are coupled by the equations of motion ([]). This coupling is of crucial 
importance in order to treat the correlation effects associated with the propagation of ele- 
mentary excitations in an interacting system. The equations of motion have also a structure 



formally identical to the one of the Bogoliubov equations for the dilute Bose gas []27| and 



to the one of the Beliaev equations for Bose superfluids [ |TSpS| . With respect to those the- 
ories the present approach makes use of a finite ranged and momentum dependent effective 
interaction. The same theory gives a reliable description of surface modes (ripplons) both 
at small and high momenta ||21|| . In the vacuum the equations of motion coincide with the 



Schrodinger equation for the free atom wave function f(z), while g(z) vanishes. 

The equations for f(z) and g(z) are solved numerically in a finite box. We work in a 
slab geometry, that is, with the liquid confined within two parallel surfaces, the slab being 
centered in the box. The box size is of the order of 100 -j- 150 A and the slab thickness 
is typically 50 -j- 100 A to simulate sufficiently well the semi-infinite medium. The space 
along z is represented by a mesh of points, with step 0.1 -r- 0.15 A. The Hartree equation 
(|j) for the ground state is solved by diffusing a trial wave function in imaginary time. After 
several iterations the procedure converges to the density profile which minimise the energy 
functional. The minimisation takes 2 -=-5 hours of CPU-time on a workstation, depending on 
the accuracy required. The same calculation provides the self-consistent mean field entering 
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the static Hamiltonian (||). Then we expand the linearised time dependent wave function 5^/ 
on the basis of eigenstates of H . By this way the equations of motion for / and g become 
a matrix equation for the coefficients of the expansion. We take a basis of 100 states or 
more to cover all the relevant part of the spectrum, so that the matrix to diagonalise is at 
least of dimension 100 x 100. The solution of the equations of motion for a given value of 
parallel wave vector q x takes about two hours of CPU-time on a workstation. The output is 
a set of discrete eigenenergies and the corresponding functions f(z) and g(z). By solving for 
different values of q x one can span the whole spectrum. An example is shown in Fig. |2], where 
we plot the energies of the excitations of a slab 80A thick in a box of 150A. One recognises 
clearly the main features already reported in Fig. [I], namely the phonon-roton dispersion 
and threshold for free atom states. The spacing in energy depends on the thickness of the 
slab and of the free space between the slab surface and the confining box. In the limit 
of infinite slab thickness the excitations in the liquid will reduce to a continuum of states, 
as already discussed for Fig. |1]. In Fig. Q one also notices lower branches below the roton 
energy A and the phonon dispersion. They correspond to surface modes and a discussion 
about them has been already reported in Ref. |21|| . The function / and g, at given q x and 
u, will take contributions from the elementary modes propagating inside and outside the 
liquid at different angles with respect to the surface. The contribution of each elementary 
mode (phonon, R + and R~ roton, atom, surface mode) is characterised by the appropriate 
dispersion law and can be identified by looking at the form of f(z) and g(z) and their Fourier 
transforms. An example is given in Fig. ^, where we show the form of f(z) and g(z) for 
one of the states of the spectrum in Fig. |2|. The amplitude and the wave vector associated 
with the free atom are easily extracted from the form of f(z) close to the boundary. The 
amplitudes and wave vectors of rotons can be calculated by looking at the Fourier transforms 
of f(z) and g(z) in the liquid. In the lower part of Fig. |3| we show the Fourier transform of 
f(z). In the limit of infinite slab thickness, the function f(q z ) would reduce to the sum of 
two (^-functions of the form f±5(q z — q z ±), corresponding to i? ± rotons. Since the Fourier 
analysis is restricted to a finite box within the slab, f{q z ) is actually the convolution of those 
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5-functions with the transform of a step function; this convolution produces the oscillations 
visible in the figure. The amplitudes f± and the wave vectors q z ± are extracted with great 
accuracy by fitting f(q z ). In the case of Fig. |3|b the best fit is indistinguishable from f(q z ). 

All these numerical results are used as input in the calculation of the flux of excitations 
in the scattering processes at the surface, as explained in the following section. 



III. FLUX OF ELEMENTARY EXCITATIONS 

In the standard theory of scattering the probability associated with a given output chan- 
nel is written as the ratio between the outgoing flux in that channel and the incoming flux. 
In the case of quantum evaporation the free surface acts as scattering region and one has to 
consider the flux of elementary excitations, i.e., rotons, phonons and atoms. The solutions of 
the equations of motion provide the fluctuations of the density and current of helium atoms 
expressed by means of functions f(z) and g(z) entering 5\F These functions are different 
inside and outside the liquid and the boundary produces the appropriate matching between 
the excitations of the bulk liquid and the free atoms states in vacuum. Now one has to relate 
the asymptotic form of these functions with the flux of incoming and outgoing excitations 
in a given scattering process. We proceed in two steps. First we introduce the concept of 
current of elementary excitations within the formalism of linearised time dependent mean 
field theory. Then we represent the scattering processes by taking linear combinations of 
stationary solutions of the TDDF equations. 

Far away from the surface and for a given energy hu the function 5^ is the sum of plane 
waves of the form 

6*(t) = £ [/^ + j;(r)e w ] = £ [fjjto*-"* + g,e- 1 ^-^} (13) 

3 ' ' j 

where the index j = a, +, — refers to the type of excitations (atoms, R + and R~ respectively) 
which contributes to <5\F Each type of excitation obeys its own dispersion law uj(q) and 
contributes to S^f with amplitudes fj and gj. 
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Well outside the liquid 5^ describes free atoms. The hole-particle part g a vanishes, while 
the atom density is given by p a — f%. The density of current is given by 

ia = ^f a =V a f a , (14) 

where v is the group velocity which, for the free atoms, coincides with ftq/m. These are the 
usual definitions of density and current density for free particles. 

Inside the liquid the same definitions can be generalised to describe the elementary 
modes of the correlated system. The RPA (or, equivalently, Bogoliubov) formalism provides 
a natural generalisation for the density of elementary excitations which, in term of the 
components / and g of the time dependent solution (|j), is given by 

,if- ./)!••) 2 uM)' 1 fj-nj- (15) 

This expression is consistent with the well known ortho-normalisation property of the RPA 
solutions: 

JdT[f;(r)f j (r)-g*(v)g j (r)]=8 ij . (16) 

The total density of elementary excitations, p ex = J2j P*j X , satisfies the equation of continuity: 

dp e 



dt 



+ V-r = • (17) 



The quantity 'f x = J2j j| x is identified with the current density of excitations, for which one 
expects the result 

jf vMi" - vji f J ■-!!]) (18) 

where Vj = V q a>(g) is the group velocity of the j-th elementary excitation. The group 
velocity should not be confused with the quantity frq/m (which is the correct expression 
only for free atoms) and actually becomes opposite to q in the R~ range of wave vectors. In 
Appendix B we show that the form flilf) for the current density is fully consistent with the 
equations for p ex derived starting from the solutions of the TDDF equations. We note also 
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that, when p ex and j ex are calculated by using the stationary solutions of the equations of 
motion, the time derivative vanishes and the equation of continuity implies simply that the 
^-component of the current of elementary excitations is equal inside and outside the liquid. 
This result is also connected with the unitarity of the scattering matrix, as we will discuss 
in the next section. 

Equations ( 15|,T7|,|18|) emphasise a remarkable feature of the TDDF solutions. In a dilute 



system they reduce to the usual expressions holding for free particles, while in the interacting 
liquid they include the correlations associated with the propagation of the collective modes. 
For instance in the long wavelength phonon regime (q — > 0) the group velocity coincides 
with the sound velocity and g — f ■ 

From the knowledge of the current of elementary excitations associated with each asymp- 
totic solution one can calculate the flux of incoming and outgoing excitations in a scattering 
process. Actually the output of our TDDF equations, as discussed in Sec. ||, are stationary 
solutions, namely real functions f(z) and g(z) which vanish at the border of the computa- 
tional box. These functions can be viewed as a combination of R + and R~ rotons in the 
liquid and atoms outside, travelling to and from the surface as shown schematically in the 
upper part of Fig. f|. Using Eq. (]TS|) one can calculate the current associated with each 
component. In order to select a single scattering process one can combine two or more 
stationary solutions at the same energy and parallel wave vector ||§. If f {1 \z), f 2 \z) and 



f( 3 '(z) are three solutions of the linearised TDDF equation at the same energy and parallel 
wave vector, then the combination 

f{z) = f 1 \z)+a 2 f {2 \z) + a^\z) (19) 

is also a solution (the same holds for g(z)). One can choose the coefficients of the combination 
in such a way that part of the incoming flux vanishes. For instance, if we want to represent 
one R + roton impinging onto the surface and reflecting as R + or R~ , or evaporating an atom, 
we have to combine the TDDF solutions in order to make the incoming fluxes of atoms and 
R~ rotons vanish. This provides two equations for the complex coefficients a 2 and 03 in 
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terms of the amplitudes and (the third coefficient a\ fixes the normalisation of the 
wave function and its value is not relevant when one calculates the ratio of fluxes). Once the 
coefficients are calculated, one is left with a solution which represents a flux (current density 
projected along z) of incoming R + rotons and a flux of evaporated atoms and reflected R + 
and R~ rotons, as in the lower part of Fig. |j. By definition, the evaporation probability is 
given by 

flux of evaporated atoms 

P+a = 7. T~. : 777 — j (20) 

flux of incoming R + 

the normal reflection probability is 

„ flux of reflected R + , . 

p ++ = a — r : — p+ ( 21 ) 

flux of incoming R + 

and the mode-change reflection is 

flux of of reflected R~ 
+ flux of incoming R + 

The explicit formulae for this scattering process, including the geometry of the slab and 
the box, are given in Appendix C. By combining the stationary solutions with different 
coefficients one can select the process in which the incident excitation is a R~ roton. In this 
case the ratios of the outgoing fluxes and the incoming R~ flux provide the probabilities 

P-a (evaporation), P (normal reflection) and P_ + (mode-change reflection). Similarly, by 

selecting the process with an incident atom one gets the probabilities P aa (reflection), P a+ 
(condensation in R + ) and P a - (condensation in R~). In region III of the spectrum in Fig. [I] 
these are the only one-to-one processes available. Generalisation to the other regions of the 
spectrum is straightforward. 

The method we have used in Ref. |B| to derive the evaporation and reflection proba- 
bilities is simpler but less general. It corresponds to the case when one of the state in the 
linear combination has vanishing amplitude outside the liquid, due to destructive interfer- 
ence. The existence of such resonant states simplifies the calculation of the probabilities P^. 
However, the more general method used in the present work is not restricted to special states 
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and provides results on a wider range of energy. Furthermore it can be easily extended to 
describe more scattering channels. 



IV. THE SCATTERING MATRIX 

The scattering matrix S characterises a scattering process through the definition ty ut = 
Sty i n . The unitarity and the time reversal symmetry (t — ► —t) of S imply 

S ] = S- 1 (23) 

and 

S* = S- 1 , (24) 

respectively. Combining these two conditions one finds that the scattering matrix elements 
must satisfy the general property 

Sij — Sji . (25) 
The scattering matrix elements are related to the probabilities P^- by 

Pij = \Sij\ 2 ■ (26) 



Equation ( p5|) then implies that the evaporation, reflection and condensation rates, intro- 
duced in the previous section, reduce to six quantities to be calculated: P +a , P_ a , P ++ , P , 

P + _ and P aa . Furthermore the unitarity condition takes also the form S*S = 1, or 

|5 , +a | 2 + |5 ++ | 2 + |5 , + _| 2 = l (27) 
\S. a \ 2 + |5__| 2 + |5_+| 2 = 1 (28) 

|5 +Q | 2 + |S_ Q | 2 + |5 aa | 2 = l. (29) 

These conditions express the conservation of flux in the three processes where the incoming 
excitation is a R + roton, a R~ roton and an atom, respectively. They reflect the fact that 
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the present approach does not include inelastic channel. The unitarity of the S matrix gives 
three additional constraints: 



+ 







(30) 







(31) 



^^n*^— a ~T O i . O — l A O. 



. 



(32) 



All these conditions for the scattering matrix elements can be used to obtain useful checks 
for the results of the numerical calculation. 

V. RESULTS 

The main results of the numerical analysis are shown in Fig. |5|, where we report the 
probabilities Pjj = |SV/| 2 as a function of energy. The data are obtained with several values 
of the parallel wave vector q x in order to span the phonon forbidden region (region III of 
the spectrum in Fig. [I]). The data at lowest energies in Fig. [| correspond to q x = 0.6 A -1 , 
the ones at highest energy to q x = 0.85 A -1 . From the knowledge of the orthogonal wave 
vector q z one extracts also the incidence angle. For R + and R~ rotons it is approximately 
15° -T- 20° and 20° -j- 30°, respectively, while for atoms is 50° -j- 90°. The range of incidence 
angles for rotons is rather narrow, so that a discussion about the angular dependence of 
the evaporation and reflection probabilities is not significant in region III of the spectrum. 
Further work is planned to explore the angular dependence including region II. The error 
bars arise from fluctuations of the numerical results. They originate mainly from the fact 
that the states in the linear combination ( |T9"D may be not enough linearly independent. 
As explained at the end of Appendix |C|, this implies an amplification of small numerical 
uncertainties in the values of amplitudes, wave vectors and other parameters entering the 
expressions for P^. 
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Looking at Fig. |5] one notices that the unitarity conditions ( p7| - p9| ) are satisfied within the 
error bars. Actually each point with error bar corresponds to several combinations of states 
at the same energy. The unitarity condition can be checked for each combination separately. 
It turns out that it is satisfied within about 5% when the states in the combination are 
sufficiently linearly independent. An example of numerical values is given in Table | for 
states at energy 10.87 K. The results of different combinations of four states are given. 
The four states correspond to symmetric and antisymmetric solutions for two slabs (62 and 
63.7 A thick). One notes that the values of the probabilities P^ extracted form different 
combinations are close each other. The unitarity conditions, which correspond to E = 1 in 
the notation of the table, is satisfied within 5% in all cases. Furthermore one has Pjj ~ Pji 
as expected from Eq. (f25p. 

The evaporation probabilities for R + and R~ rotons behave quite differently. The one 
for R + rotons grows from to 1 by increasing the energy. A similar behavior has been 



recently observed in experiments at normal incidence in the same range of energy pO 
Conversely the R~ rotons are less effective in evaporating atoms. The ratio between the two 
probabilities has to be 1 at the energy A of the roton minimum where the two excitations 
become identical. As soon as the energy grows the ratio P +a /P_ a increases and becomes 
large. Results for P +a /P_ a are shown in Fig. ||. 

At low energies, approaching A, the mode-change is dominant in the reflection processes 
for rotons. This is also consistent with symmetry arguments (see discussions in Ref. [pl[]). 
At higher energies the mode-change becomes less probable. An incident R~ rotons has low 
probability of evaporating atoms, as discussed above, and it reflects as R + at low energy and 
as R~ at high energy. The behavior of an incident R + roton at high energy is completely 
different: the normal reflection as R + is found to be practically zero everywhere. R + rotons 
prefer to evaporate atoms or to reflect as R~ , the sum of the two probabilities being almost 
1. This fact can be explained qualitatively by considering the change of momentum in a 
normal mode reflection. This should be twice hq z , where q z in this part of the spectrum is of 
the order of 2 A -1 . But it is difficult to imagine a mechanism for absorbing a momentum of 



the order of 4 A -1 at the surface, the latter being quite smooth (the Fourier transform of the 
density profile has vanishing components above 2 A -1 ). On the contrary, the momentum 
transfer for reflections including R~~ rotons is at least a factor two smaller and it is expected 
to be more probable. 

The fact that P ++ is almost zero everywhere has an interesting consequence. If we assume 
S ++ = and use the unitarity conditions (|27|-|32j) , we can express all the probabilities through 
a single parameter, for instance P +a . Simple algebra gives (see Appendix D): 

F+- = 1 - P +a (33) 

P- = P+a (34) 
Paa = Pi- = (1 " P+a) 2 (35) 
P_ a = P-+P+ a = (1 " P+a)P+a (36) 

Note that the maximum value of P- a , from the last equation, is obtained for P +a = 1/2 and 
thus the following inequality holds: 

P-a < \ ■ (37) 

All the above relations, which follow from the assumption S ++ = 0, are well satisfied by the 
numerical results in Fig. |^. 

Finally, the results for the reflection probability for incident atoms are more intriguing. 
One knows from the experiments that the reflection probability is less than 1% apart from 
very small values of q z where the long range Van der Waals interaction dominates. Our 
results for the reflection probability are consistent with the experimental results only at 
high energy, where P aa drops to zero. At lower energy the roton channels become less and 
less active in the condensation process and the atom is reflected with large probability. This 
behavior is mainly the consequence of the fact that in the range of angles considered in the 
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present work phonons cannot be excited. As soon as phonons are allowed (region I and 
II in Fig. [I]) they are expected to provide the natural channel for the condensation of the 
incoming atoms. It has been also suggested that, as soon as the atom has enough energy to 
feel the external part of the density profile, it dissipates its energy by producing many low 
energy ripplons. This mechanism is not accounted for in our one-to-one theory. 



VI. CONCLUSIONS AND PERSPECTIVES 

In this paper we have investigated the evaporation and condensation processes taking 
place at the free surface of superfluid helium. The main achievements of the present work 
concern both the formal development of scattering theory, in the framework of linearised 
time dependent density functional formalism, and the numerical analysis, which provides 
first systematic results for the evaporation, condensation and reflection rates. The most 
relevant steps and achievements of our approach are here summarised: 

• i) The equations of motion of linearised time dependent density functional theory have 
been explicitly solved in the presence of the free surface using the Orsay-Trento phe- 
nomenological density functional. By taking linear combinations of different solutions 
carrying the same energy and momentum parallel to the surface, we have selected the 
relevant scattering processes (incoming atom reflected into an atom or condensed into 
a R + or R~ roton, etc.). 

• ii) We have identified, in the framework of the time dependent density functional 
scheme, the current of elementary excitations which obeys the equation of continuity 
( |17D and allows one to evaluate the rates of evaporation and reflection. This current 
is determined by the group velocity and by the density of elementary excitations for 



which the familiar RPA expression ( |T5| ) holds. The resulting formalism has been shown 
to be consistent with the general properties of the scattering matrix, following from 
unitarity and time reversal invariance. 
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• iii) The numerical results, restricted in the present work to the region where only 
condensation of atoms into rotons is possible (region III of Fig. |1|), reveal that rotons 
R + are systematically more effective than rotons R~ in the evaporation mechanism. An 
interesting feature emerging from our analysis, is that the P ++ scattering probability 
is extremely small at all energies. This behavior, with the help of the formal properties 
of the scattering matrix, allows one to express all the evaporation and reflection rates 
in terms of a single parameter [see Eqs. ( |33l]36|) 1. 

• iv) A major discrepancy between our predictions and the experimental data concerns 
the atomic reflectivity which is predicted to be significantly large when the energy 
approaches the roton minimum. This behavior is easily understood since our theory 
accounts only for one-to-one processes where condensation into ripplons is excluded by 
the conservation laws. Due to the absence of phonon excitations, which are excluded in 
region III, no active excitations are allowed in the condensation process when the en- 
ergy approaches the roton minimum with the consequent full reflection of the incident 
atom. 

Natural developments of the present work will include: 

• The study of the scattering process in different ranges of energy and incident angles, 
where phonons can participate in the scattering process. 

• The study of multi-ripplon excitations whose effects can be taken into account through 
the use of an optical-type potential in the time dependent equations. 

Work in these directions is in progress. 
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APPENDIX A: ORSAY-TRENTO DENSITY FUNCTIONAL 



We use the phenomenological density functional of Ref . p(J . It has the form: 



E = £ (fcm) [p, v] + [p] + E^ [p, v] , ( Al) 

where p and v are the density and velocity of the atoms, respectively. The first term is the 
kinetic energy of the non interacting bosons, 

h 2 



E^ km \pM = [dr — |V^(r)| 2 
J 2m 

= /rfr{|l(VVp) 2 + ^p(r)|v(r)| 2 } ; (A2) 

the correlation energy 

E (c) 

is given by 

E (c) [p] = Jdv{\jdv> pm^-r'Mv') + |p(r)(p r ) 2 + |p(r)(p r ) 3 
-a. I dr' F{\v - r'|) ( 1 - Vp(r) • Vp(r') ( 1 - } ; 



A 2 

4m 



(A3) 



finally, the backflow energy 

E (bf) 

is 



EM\p,v] = -j J jdvdv'Vj{\v-v'\) p{v)p{v>) [v(r)-v(r')] 2 . (A4) 

The first term in the kinetic energy, which depends on gradient of the density, is a quantum 
pressure; it corresponds to the zero temperature kinetic energy of non-interacting bosons of 
mass m. The two-body interaction Vi in the correlation energy E^ is the Lennard- Jones 
interatomic potential, with the standard parameters a = 2.556 A and e = 10.22 K, screened 
at short distance (V = for r < h, with h = 2.1903A). The two terms with the parameters 
c 2 = -2.411857 x 10 4 K A 6 and c 3 = 1.858496 x 10 6 K A 9 account phenomenologically 
for short range correlations between atoms. The weighted density p is the average of p(r) 
over a sphere of radius h. Those terms are very similar to the functional of Ref. p4j. The 
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last term in E^ c \ depending on the gradient of the density in different points, has been 
added in order to improve the description of the static response function in the roton region. 
The function F is a simple Gaussian, F(r) = 7r~ 3//2 £ -3 exp(— r 2 /£ 2 ) with I — 1 A, while 
a s = 54.31 A 3 and po s = 0.04 A~ 3 . The energy E^^ contains an effective current-current 
interaction accounting for backflow-like correlations. In Ref. |2(J the simple parametrisation 

Vj{r) = (7xi + 7i 2 r 2 ) expf-a^r 2 ) + (721 + 7 22 r 2 ) exp(-a 2 r 2 ) , (A5) 

was chosen in order to reproduce the phonon-roton dispersion in bulk liquid. The parameters 
are given in Table [TT]. 

APPENDIX B: EQUATION OF CONTINUITY AND CURRENT OF 
ELEMENTARY EXCITATIONS 

In this appendix we show that definition ( |18|) for the current of elementary excitations 
is consistent with equation (||) holding in time dependent density functional theory. In 
order to make the demonstration more transparent we use a simplified version of the energy 
functional of the form 

H = |h W( ri )| 2 + J dr 2 V(\r 1 - r 2 |)p(n)p(r 2 ) , (Bl) 

where V is a generic two-body interaction. The generalisation to the more complex Orsay- 
Trento functional is conceptually straightforward but rather tedious in practice. 
Taking variations with respect ty* in Eq. (f|) one gets 

*(r l5 t)=0. (B2) 

Here we are interested in the solutions well inside the liquid, where the ground state density 
p is a constant and the above equation yields 

li = Pol drV(r) . (B3) 
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As in Eq. (|]) we take a linear expansion of \I/ to describe the excited state. Let us include 
the time dependence in the definitions of / and g in the form 

f(r,t) = fe 1 ^-^ ; g(r, t) = ge 1 ^'^ , (B4) 

so that 

6*{r,t) = f{T,t)+g*(T,t). (B5) 



Inserting this expression in Eq. (|B2|) and isolating the terms oscillating with positive and 
negative frequencies, one gets two coupled equations for / and g: 



ih 



0/(1) & 
dt 



— V 2 /(l)+Po Jdv 2 V(r 12 )[f(2) + g(2)] = (B6) 



ih 



dg(l) h 2 
dt 2m 



V 2 g(l) + po J dv 2 V(r 12 )[f{2) + ^(2)] = (B7) 



where we have used the short notation (1) = (ri,i) and (2) = (r 2 ,t). Now, take /(l) times 
the complex conjugate of Eq. ( |B6|) and subtract /*(1) times Eq. (|B6|) itself. Similarly take 
g(l) times the complex conjugate of Eq. ([B7]) and substract g*(l) times Eq. (|B7p . Summing 
up all the terms, one gets 

*^(i/(i)i 2 - i^(i)i 2 ) = ^ [/(i)v 2 r (i) - r (i)v 2 /(i) + <?(i)vy (i) - g*(i)v 2 g(^ 



- po J dr 2 V(r 12 ){f(l)f*(2) - /* (l)/(2) + f(l)g*(2) - g*(l)f(2) 
+ g{l)f*{2) - r 0)9(2) + g(l)g*(2) - ^(1)^(2)1 . (B8) 



This equation has the form of an equation of continuity if the quantity 

ff* = \f(l)\* - \g(l)\* = f* - <? (B9) 

is interpreted as the density of elementary excitations and if the right hand side is identified 
with —ifiV ■ 3 ex . Now we prove that the quantity j 6X has indeed the form 'f x = vp ex , where 
v is the group velocity of the excitation. To this purpose we rewrite properly the various 
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terms in Eq. (|B8|). We begin with the one coming from the kinetic energy. All terms with 
the laplacian can be written in the form —ihV ■ [j ex ]kin, with 

[TV = [/(l)Vf(l) - /*(1)V/(1) + g(l)Vg*(l) - g*(l)Vg(l)] . (BIO) 

We remember that /(l) and g(l) are plane waves, as in Eq. (|B4|). Thus the expression for 
IT^fcm, in terms of the real amplitudes / and g, becomes 

[rw=^V + <? 2 ), (bii) 

To manage the terms containing the interaction is less trivial. Let us use the Fourier de- 
composition 

/(r,t) = Jdqxf{qx,t)e^ (B12) 

and take the limit /(qi) = /5(qi — q) at the end. We can rewrite the first two terms with 
the interaction in Eq. (|B8|) as 



/dr 2 V(r 12 )|/(l)r(2)-r(l)/(2)} = J J dq^ e^~^ 

x /(qi,t)r(q 2 ,t)(V(q 2 ) - V^)) (B13) 

where V(q) is the Fourier transform of the interaction. By doing the same with the terms 
containing g, one gets 

- iftV ■ [f x U = -po J /rfqiq 2 e i(qi - q2)ri (^(q 2 ) - V(qx)) [f(qi, t)f* (q* t) 

+ <?(qi, t)/(q 2 , *) + /(qi, %*(q 2 , *) + <?(qi, *)/* (<J2, *)] • (B14) 

Now, by considering /(q) and g(q) as delta functions, one gets 

LFW = Po(/ + <?)(/* + 9*)V q V(q) • (B15) 
and summing [j ex ] p0 i to [j e:r ]fcm one finally obtains the result 

[ (f + g 2 ) + Po(f + g)(F + gT)V 9 v(q) (Bie) 
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m 



for the total current. At this point it is convenient to rewrite Eq. ( P16|) in terms of the 
dispersion law u(q) coming from the solutions of the equations of motion ( |B~6| - |B~7| ) which, 
for plane waves of the form (fB3j) and after some algebra, yield 

hq 2 



f + 9 



2muj 



(f-g) 



(B17) 



with 



f + 9 2 



2m 
Tiq 2 uo 



10 



'hq^ 
2m 



(f 



2 <f: 



Inserting these results in Eq. ( |B16|) one finally finds 



(B18) 



(B19) 



j ex = v(/ 2 - g 2 ) = vp 



(B20) 



where v coincides with the usual definition of the group velocity: 



v = V q u(q) 



(B21) 



This is the desired definition of the current density, in terms of which Eq. 
of the equation of continuity (|TT|): 

dp ex 



dt 



+ V-j e 



. 



has the form 



(B22) 



APPENDIX C: DETERMINATION OF THE SCATTERING RATES 

We give an example of derivation of the probabilities Py in the case of R + rotons coming 
to the surface and producing evaporated atoms and reflected R + and R~ rotons. The 
geometry we use is the one in Fig. [j]: the surface is at z = 0, the center of the slab at z = L 
and the left border of the box at z — —b. The solutions of our TDDF equations are real 
functions which vanish at z — —b. They can be either symmetric and antisymmetric with 
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respect to the center of the slab (z = L). If b and L are sufficiently large, the functions f(z) 
and g(z), characterising the changes 5^ as in Eq. (^), close to the boundary and near the 
origin are simple sinusoidal functions. In particular, close to z = — b one has a free atom 
state of the form 



f(z) ~ /„ sm[q a (z + b)] = ^. e iqab e iqaZ - ^ e ~ iqal 'e~ iqaZ 



(CI) 



and g = 0. Close to z = L, inside the liquid, f(z) is the sum of two oscillating functions of 
the form 



f(z) ~ /± sm[q ± (z - L)\ = f -± e ~^ L e^ z - f -± e iq ± L e~ iq ± z 



(C2) 



in the antisymmetric case (a cosinus in the symmetric case) and the same for g. Here, 
for simplicity, the symbol q, for wave vectors, denotes the ^-component only. The Fourier 
analysis of the solutions of the TDDF equations provides the amplitudes f a and f± for fixed 
values of L and b. Changing L and b one can find different solutions at the same energy and 
parallel wave vector. Let us choose two symmetric states and an antisymmetric state and 
combine them in the form 



f = f {1 \z) + a 2 fW(z)+a 3 fW(z) 



(2), 



f(3), 



(C3) 



where the solution corresponds to the parameters L a and b a . The atomic part, well 
outside the slab, becomes 



fa 



2i 



(2) iq a b 2 



-tq a z 



2i 



fV) e -«J* + a 2 fWe-^ + a 3 f?>e 



(2) -iq a b 2 



-(3)„-ig<ii>3 



(C4) 



the R + roton part, close to z = L, is 



A 



-if^e iq+bl + a 2 f { + ] e iq+h2 + a 3 f^e iq + b3 



-iq+z 



+ 



if(X) e -iq+W + a2 f^ e -il + b2 + a ^fVV e 



(3) -iq + b 3 



(C5) 



and the R roton part is 
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/- 



iq-Z 



-tq—Z 



+ 



-if { y e iq - bl + a 2 f i L ) e iq - b2 + a 3 f^e iq - b3 
if ( L ) e- iq - bl + a 2 f i - ) e- iq - b2 + a 3 f { * ] e~ iq - bz 



(C6) 



The quantities in brackets represent the amplitudes of incoming and outgoing excitations 
and one can choose the coefficients a 2 and a 3 in order to make two of them vanish. As said 
at the beginning, we show here the case of an incident roton R + . So, we must have zero 
amplitude for incoming atoms and Rr rotons (see Fig. ^). Keeping in mind that R~ rotons 
travel with negative group velocity, one obtains the following algebraic equations for a 2 and 
a 3 : 



if m e- iq - bl + a 2 f ( J'e- iq ' b2 + aj^'e 



( 2 ) tD -iq-b 2 



c(3) -iq- 63 



. 



(C7) 



Since the ratio between the amplitude of the particle-hole / and the hole-particle g compo- 
nents depends only on the energy and not on the particular state in the combination, the 
above equations ensure that the incoming flux of atoms and R~ rotons vanish. One has 



to recall the definition (|T8|) of the current density and use the values of a 2 and 03 coming 
from Eqs. ( |CT| - |C8| ) in order to calculate the flux of the incoming rotons and the ones of the 
outgoing excitations. The scattering rates come out to be: 

-iA { + ] e iq + bl + a 2 A { l ] e iq + b2 + a 3 A^e iq + ba 2 



iA (l) e - iq _ bl + ^m^iq-bz + a A {3) e -iq_b 3 



iA i l ) e- iq + bl + a 2 A ( + ] e~ iq + b2 + asA^e-^^ 



A(i) e -iq a bi + a 2 Af) e - iqab2 + a 3 A^e- iqab3 



< e -iq+h + a% j££ > e -iq+b 2 + a 3 A% } e- iq + b * 



I (2) 



(3), 



(C9) 



(CIO) 



(Cll) 



where Aj = sign(fj)yvj(ff — gj). One can also derive similar expressions using two anti- 
symmetric states and one symmetric state. The other scattering rates can be evaluated by 
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selecting different processes, i.e., by solving different equations for the coefficients a 2 and 
a 3 . The states entering the combinations have to be linearly independent. This requires 
a certain care in selecting proper values of L and b. If the states are not sufficiently inde- 
pendent, it happens that the denominators in the expressions for Pij tend to be small and 
consequently small fluctuations in the numerical results produce large uncertainties in the 
final results. 

APPENDIX D: UNITARITY CONDITION IN THE ABSENCE OF R + NORMAL 

REFLECTION 



In this appendix we will derive expressions ( p3| - pq ) for the scattering matrix elements in 
the phonon forbidden region (region III in the spectrum of Fig. [I]) which are valid if S++=0. 
In this hypothesis one can rewrite Eqs. Q27|), Q5DD and respectively, as 



\S +a \ 2 + \S+-\ 2 = l (Dl) 

S* aa S +a + S* a _S + - = (D2) 

S* +a S^ a + Si_S— = . (D3) 
The first equation coincides with the result (P3). From the other two one gets 



S1LS- = S* + 2 a S aa (D4) 



or 



1 5 -| — \^\S — | 2 — |5' +a | 4 |S' aa | 2 . (D5) 
Subtracting Eq. ( Effi) from Eq. (|28| ) one gets 

\S-\ — | 2 — I'S'+aP = I'S'aal 2 — l^* | 2 . (D6) 



Excluding S aa from Eqs. ( |D5|) and (|D6| ) one finds 
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|<S+-| 2 — I ^S+a. | 2 



\S- 



\S+a\ 



i)\s-\ 



(D7) 



or, using ( pi|) : 

|S__| 2 



\S+a 






2 — \S+a 


2 ) 


\S+- 


4 _ 


' \S+a 


4 



\S + a 


■I 


\S+- 


2 + \S + a\ 2 



\3. 



+a\ 



(D8) 



Thus we have obtained expression (B3). Expression can be derived in analogous way. 



Finally, by substituting flDl|) and ( p8|) into (|28|) , one finds 



|5_ a | 2 = l-|5 + _| 2 -|5__| 2 

= |*S'-f-ot| 2 — |*S+a| 



(D9) 



which corresponds to the last result given in Eq. (pq). 
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FIGURES 

FIG. 1. Schematic picture for the spectrum of elementary excitations in presence of a free 
surface. Solid line: phonon-roton dispersion for excitations propagating in the direction parallel 
to the surface (q z = 0). Dashed line: threshold for atom evaporation. Excitations propagating at 
different angles {q z 7^ 0) fill regions I- VI, as explained in the text. 

FIG. 2. Spectrum of elementary excitations for a slab 80A thick in a box 150A wide, as obtained 
solving numerically the TDDF equations. Two surface modes (ripplons) are also visible below the 
phonon-roton dispersion. For simplicity, only states symmetric with respect to the center of the 
slab are shown. 

FIG. 3. Example of functions f(z), in arbitrary units, and g(z), in the same units but multiplied 
by 5, for a state at energy 10.2 K and parallel wave vector q x = 0.7A -1 in the same slab as in Fig. ^. 
The density profile of the slab is shown as a dot-dashed line. In part b, the Fourier transform of 
f(z) is shown; it corresponds to the convolution of two delta functions, at the R + and R~ roton 
wave vectors, with the transform of the step function used to restrict the Fourier analysis within 
the slab. The function g(q z ), not shown, has a similar structure, but much smaller amplitude. 

FIG. 4. Schematic picture for scattering processes at the free surface. In the lower part a par- 
ticular process is selected: an incoming R + roton can produce atom evaporation, normal reflection 
or mode-change reflection. 

FIG. 5. Results for probabilities Pij as a function of energy, in the phonon forbidden region. 

FIG. 6. Ratio of the evaporation probabilities for R + and R~~ rotons as a function of energy. 
Open circles: q x = 0.6A -1 ; solid circles: q x = 0.65A -1 ; open squares: q x = 0.7A- 1 ; solid squares: 
q x = 0.75A -1 ; triangles: q x = 0.8A -1 . Dashed line: limiting value P +a /P- a = 1 for energy 
approaching A = 8.6K. 

FIG. 7. Schematic picture for the geometry used in Appendix C. 
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TABLES 

TABLE I. Probabilities Pjj for three different scattering processes (incident atom, R + and R~ , 
respectively) described by linear combinations of four states at the same energy, 10.87 K, and the 
same parallel wave vector, 0.7A -1 . The states 1-4 are symmetric (2 and 4) and antisymmetric (1 
and 3) solutions for L slab = 62A (1 and 2) and 63.7A(3 and 4), and for L box = 126A (3), 134A (4) 
and 135A (1 and 2). The quantity S is the sum of the probabilities for each process. 





atom — ► atom, R + , R 


R + -> atom, R+, R~ 


R — ► atom, R + , R 


comb. 


Paa P(X+ Pa- £ 


P++ P+- P+a £ 


P— P + P-a £ 


1-2-3 
1-2-4 
1-3-4 
3-2-4 


0.148 0.654 0.224 1.026 
0.146 0.645 0.224 1.015 
0.148 0.650 0.223 1.021 
0.146 0.645 0.224 1.015 


0.0004 0.357 0.603 0.960 
0.0003 0.358 0.612 0.970 
0.0004 0.358 0.605 0.963 
0.0004 0.358 0.612 0.970 


0.404 0.387 0.224 1.015 
0.403 0.385 0.228 1.016 
0.403 0.387 0.225 1.015 
0.404 0.384 0.228 1.016 



TABLE II. Values of the parameters used in Vj(r), see Eq. ([A 4}). 



711 


721 


712 


722 


a\ 


«2 


-19.7544 


-0.2395 


12.5616 A" 2 


0.0312 A" 2 


1.023 A" 2 


0.14912 A" 2 
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